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ABSTRACT 

The Swift/BAT 9-month catalogue of active galactic nuclei (AGN) provides an unbiased cen- 
sus of local supermassive black hole accretion, and probes to all but the highest levels of 
absorption in AGN. We explore a method for characterising the bolometric output of both 
obscured and unobscured AGN by combining the hard X-ray data from the SwiftfBAI instru- 
ment (14-195keV) with the reprocessed IR emission as seen with the Infrared Astronomical 
Satellite (IRAS) all-sky surveys. This approach bypasses the complex modifications to the 
SED introduced by absorption in the optical, UV and 0.1-10 keV regimes and provides a 
long-term, average picture of the bolometric output of these sources. We broadly follow the 
approach of Pozzi et al. for calculating the bolometric luminosities by adding nuclear IR and 
hard X-ray luminosities, and consider different approaches for removing non-nuclear contam- 
ination in the large- aperture IRAS fluxes. Using mass estimates from the black hole mass-host 
galaxy bulge luminosity relation, we present the Eddington ratios X^dd and 2-10 keV bolo- 
metric corrections for a subsample of 63 AGN (35 obscured and 28 unobscured) from the 
Swift/BAT catalogue, and confirm previous indications of a low Eddington ratio distribution 
for both samples. Importantly, we find a tendency for low bolometric corrections (typically 
10-30) for the obscured AGN in the sample (with a possible rise from ~ 15 for AEdd < 0.03 
to ~ 32 above this), providing a hitherto unseen window onto accretion processes in this class 
of AGN. This finding is of key importance in calculating the expected local black hole mass 
density from the X-ray background since it is composed of emission from a significant popu- 
lation of such obscured AGN. Analogous studies with high resolution IR data and a range of 
alternative models for the torus emission will form useful future extensions to this work. 

Key words: black hole physics - galaxies: active - galaxies: Seyfert 



1 INTRODUCTION 

Accretion onto supermassive black holes (with masses of 10^ — 
10^ solar masses) is responsible for the power output of active 
galactic nuclei (AGN). Thermal emission from an accretion disc 
emerges in t he optical-UV regime and manifests as the 'Big 
Blue Bump ' (IShieldslll978h . and inverse Compton scattering of 
UV disc photons b y a corona above the disc is responsible for 
the X-ray emission dZdziarski et al.|[l99Qb . This optical-to-X-ray 
emission constitutes the bulk of the intrinsic accretion emission 
in AGN. In the standard paradigm for AGN, a dusty torus sur- 
rounding the accretion disc partially absorbs this emission and re- 
radiates it in the infrared ((Pie r & Krolik 1992), manifesting as an 
IR bump (lElvis etal.|[l994r If the AGN is viewed through the 
torus, the optical-to-X-ray emission may be significantly absorbed 
along the line-of- sight, presenting difficulties in determining the 



true, bolometric accretion lumin osit y Lboi from the central en - 
gine. I Vasudevan & FabianI J2007h and lVasudevan & FabianI (l2009l) 
present Lboi (along with bolometric conversion ratios and accre- 
tion rates) determined from optical-to-X-ray spectral energy dis- 
tributions (SEDs) for AGN, but both studies highlight the difficul- 
ties associated with determining Lboi when spectral complexity as- 
sociated with absorption is present. In optical surveys of AGN, it 
is possible for objects to be missed altogether because such ab- 
sorption can reduce the flux severely enough (reddening) to be 
below the detection limit. Although X-ray su rveys are more ef - 
fective at providing a census of AGN activity (lMushotzkvll2004h . 
significant columns of dusty gas (more than ^ lO^^cm"^) will 
heavily reduce the X-ray (0.1-10 keV) flux and cause similar prob- 
lems. Heavily Compton-thick objects (those with column densities 
Nu ^ lO^^cm"^) in particular can have X-ray fluxes reduced by 
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many orders of magnit ude from their intrinsic flux in this energy 
band dlkeda et alj[2^ . 

The Burst Alert Telescope (BAT) on board the Swift satel- 
lite is extremely useful at addressing these considerations to pro- 
vide a more complete census of AGN activity. The very hard X-ray 
bandpass of this instrument (14-195 keV) misses absorption signa- 
tures for moderate-to-high levels of absorption. The 9-month cat- 
alogue of BAT-det ected AGN (hereafter the Swiftl'QNI catalogue, 
Tueller et al ] 12008) therefore provides an unprecedented level of 
completeness when surveying the AGN population, since it con- 
tains all but the most heavily absorbed objects. The X-ray proper- 
ties of the 153 AGN in the catalogue are presented in I Winter et al.l 
(l2009h , using data from a variety of X-ray missions to deter- 
mine their intrinsic luminosities, spectral shapes (including mea- 
sures of the spectral co mplexity) and absorbing column densities 
(A/'h). Fvasudevan et al.1 (2009) present simultaneous optical-to-X- 
ray SEDs from Swiff s X-ray telescope (XRT) and UV-optical tele- 
scope (UVOT) and calculate the total accretion output Lboi and 
associated quantities for a subsample of 26 local, low-absorption, 
low- spectral-complexity AGN from the Swift/BAI catalogue. For 
their subsample of objects, the aforementioned problems associ- 
ated with absorption are minimised or straightforwardly corrected 
for. Their work provides hints that the local AGN population is 
dominated by AGN in which X-ray emission processes contribute 
significantly to the bolometric emission, manifesting as low bolo- 
metric corrections A>:2-iokeV (^2-iokeV = i^boi/^2-iokeV ^ 
10 — 30 for 2-10 keV luminosity L2-iokev) and also suggests 
that accretion rates are low (Eddington ratios AEdd ^0.1, where 
AEdd = Lhoi/L-Edd for Eddington luminosity LEdd = 1.3 x 
10^^{Mbu/Mq) for a black hole of mass Mbh). While such an 
approach is of particular utility for determining the total power out- 
put of unobscured AGN, an alternative approach is needed to ad- 
dress this question more generally for AGN of all absorption levels. 
This study explores a method for crossing the 'absorption barrier', 
and aims to characterise the bolometric output of local AGN in the 
Swift/BAT 9-month catalogue sources across the range of absorp- 
tion properties probed by the catalogue. 

The hard X-ray data from the BAT instrument gathered over 
many months are available for the Swift/BAI catalogue sources and 
are presented online. These data provide a useful starting point for 
our absorption-unbiased study on the bolometric AGN power out- 
put. However, instead of using optical-UV emission to constrain 
the thermal emission from the accretion disc, we use the repro- 
cessed infrared (IR) emission to estimate the non-X-ray component 
of the bolometric luminosity Lboi- We take the re-processed in- 
frared (IR) emission to be a severely averaged proxy for the intrin- 
sic disc emission, since the timescale for transferral of energy from 
the high energy flux at the inner edge of the torus to emission as IR 
at i ts outer edge is of the order of several years in standard mod- 
els ( iBarvaini whereas optical, UV and X-ray variability in 
AGN is known to occur on far shorter timescales (a few minutes to 
hours). The averaged accretion output represented by the IR can be 
used in conjunction with the long-term hard X-ray observations to 
estimate Lboi- We use the archival mid-to-far infrared observations 
from the Infrared Astronomical Satellite (IRAS) all- sky catalogues 
to measure this reprocessed output; mid-infrared (MIR) observa- 
tions in particular are well suited to our purposes since they capture 
the emission in the range where reprocessed AGN emission is ex- 
pected to peak (^ 12 /am) and the all-sky surveys readily available 
in the archives are ideal for identifying IR counterparts to X-ray 
sources in the all-sky Swift/BAI catalogue. Our meth od is a low- 
redshift, absorption- 'neutral' extension of the work of IPozzi et al.l 



(2007) on high-redshift luminous obscured quasars, with the key 
difference being the use of hard X-ray BAT data instead of 2-10 
keV X-ray data. 

When considering the effects of absorption in AGN, it is im- 
portant to note that the absorbing columns measured in differ- 
ent wavebands do not always match. Under the standard unified 
scheme for AGN discussed above, one would expect an AGN with 
significant X-ray gas absorption to exhibit significant dust absorp- 
tion in the optical-UV (if one assumes a relatively constant dust- 
to-gas ratio across the AGN population). However, there are many 
classes o f object which do not obe y this simple picture, as dis- 
cussed bv lMaiolino & Risalitil ( l2007b . However, statistically speak- 
ing. X-ray and optical absorptions are broadly correlated, and sig- 
nificant mismatch occurs in only about 10-20 per cent of the AGN 
populations studied in various surveys (e.g. iGarcet et"aD l2007l . 
iMateos et "^1 120051 . [Silverman et al.|[2005h . In any case, these con- 
siderations do not not significantly affect the approach used here, 
but we return to this issue when calculating the reprocessed IR 
emission. In this paper, whenever the labels 'absorbed/obscured' or 
'unabsorbed/unobscured' are used, they consistently refer to X-ray 
absorption. 

The use of the IR emission as an indicatio n of AGN ac - 
tivity is supported by numerous works. Mushotzkv et al ] (l2008h 
find a strong correlation between the hard X-ray BAT luminosity 
and the Two-Mi c ron A ll- Sky Survey (2MASS) J-band luminosity. 
iMcKernan et al.l (l2009h use analyses of X-ray observations from 
the literature and IRAS data to calculate X-ray-to-IR luminosity 
ratios, and identify from these that the parameters specifying the 
physical AGN environment su ch as the torus geometry must span 
a narrow range. The study of iGandhi et all (l2009h discusses the 
correlation between X-ray (2-10 keV) emission and IR emission 
in detail for a sample of Seyfert nuclei, again using X-ray analy- 
ses from the literature and new IR data from the Very Large Tele- 
scope's Imager and Spectrometer for mid-Infrared (VISIR/VLT), 
taken specifically for addressing the issue of nuclear emission in 
local Seyferts. Their data are the best estimates of the nuclear (non- 
stellar) IR flux in AGN to date, and show a strong correlation be- 
tween intrinsic, uncontaminated nuclear IR monochromatic lumi- 
nosity (at 12.3 fim) and 2-10 keV X-ray luminosity. Their work 
reinforces the idea that the uncontaminated MIR continuum is an 
accurate proxy for the intrinsic power of AGN, and reinforces and 
improves constraints on previous de terminations of the correlation 
jHorst et al.ir2008l . iHorst et al.ll2006h . They also higWight the use- 
fulness of extending this approach to a very hard X-ray selected 
sample such as the one used here. Their correlation also provides 
one way to estimate non-nuclear contamination in IR observations 
taken with larger apertures, as is the case with the IRAS data used 
here. The other method employed for correcting for non-nuclear 
contamination in this study is fitting of host galaxy and nuclear 
SED templates to the mid-to-far IR SEDs. 

In this paper, we use the IRAS and BAT data to generate 
SEDs from which the reprocessed IR and intrinsic X-ray emission 
can be determined, for calculation of Lboi- We calibrate this ap- 
proach a gainst integration of the optical-to-X-ray SEDs as pre- 
sented in I Vasudevan et al.l (l2009h , using the sources overlapping 
between the two studies for this purpose. We note that although 
our IRAS and BAT data are not contemporaneous, they both rep- 
resent an averaged, long-term picture of accret ion whereas the 
simult aneous optical-to-X-ray data presented in IVasudevan et al.l 
(I2OO9I) are a 'snapshot' of the accretion emission at a particu- 
lar time;_rha^_a22roa^^ advantages (as discussed 
in I Vasudevan & Fabianll2009b . but the average, long-term picture 
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should be found to agree with 'snapshots' when considering a 
reasonably large sample of objects, since the effects of variabil- 
ity should be averaged out over many sources. Certainly, 'quasi- 
simultaneous' data between hard X-ray and IR are preferable since 
ideally we want the IR to capture a reprocessed version of the same 
long-term state seen at higher energies, but the IR and hard X-ray 
bands can afford to be separated by much longer time periods to 
qualify as 'quasi- simultaneous' due to the large timescales involved 
in reprocessing. We present the bolometric luminosities, Edding- 
ton ratios and bolometric corrections determined for all the local 
{z < 0.1) objects in the Swift/BAI catalogue with good quality 
IRAS and BAT data, using black hole mass estimates calculated 
from K-band host bul ge luminosity estimates via the method of 
IVasudevanetaD(l2009h to determine accretion rates (apart from for 
Cyg A, where a known dynamical mass estimate is used). 

This study can of fer perspectives on unification schemes for 
AGN dAntonuccilllOOi) . as it provides an idea of which type of ac- 
cretion processes dominate in AGN (as parameterised by the X-ray 
bolometric correction) over a range of absorbing column densities 
and hence orientations. A flat Universe (Qm + = 1) with a 
cosmology Ho = 71km s~^Mpc~^ and Qm = 0.27 has been 
assumed throughout this work. 



2 SAMPLE SELECTION 

We firstly apply a redshift cut ( z < 0.1) on the 145 sources with 
spectral properties presented in IWinter et iD (l2009l) and exclude 
sources which would exhibit prominent jet c ontributions to thei r 
luminosity such as BL Lacs or blazars (see IWinter et al.|[2009h , 
yielding 116 potential objects. We then exclude the two merging 
galaxies NGC 6921 and MCG -h04-48-002 since AGN at such close 
proximity are very likely to be confused in the low-resolution IRAS 
images. We also exclude the galaxy NGC 6814 which is known to 
have a prominent cataclysmic variable in the foreground, yielding 
113 potential objects. 

In assembling IRAS data for this sample, we take care to en- 
sure confident positional matches between the IRAS catalogues and 
the NED positions. We preferentially use the IRAS Point Source 
Catalogue (PSC) as our source for IR photometry; if data was not 
available in this catalogue we then turned to the IRAS Faint Source 
Catalogue (FSC), which also catalogues point sources. There was 
some overlap between the PSC and FSC; in many cases the FSC 
provided greater wavelength coverage and data quality than the 
PSC. The 12, 25, 60 and 100 /im photometry for these catalogues 
are available online along with useful diagnostic and data quality 
information. The positional accuracy of the IRAS photometry is de- 
pendent on a number of factors. We first of all restrict our searches 
to within 180 arcsec of the NED positions of our AGN (bearing 
in mind the PSF FWHM of 1-5 arcmin for IRAS), yielding 
70 matches in the PSC and 73 matches in the FSC (with signif- 
icant overlap between the two). We then employ the information 
on the positional error ellipses for the IRAS photometry provided 
for each object, to gauge the goodness of the positional match. The 
semimajor and semiminor axes (fimaj, cTmin) for each error ellipse 
are provided (in arcsec), and we determine an error circle of radius 
c^pos = (cTmaj CTmin) with thc samc area as the ellipse defined 
by these axes. We require the IRAS photometry to have a position 
within ScTpos of the NED position, taking care to ensure that only 
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one IRAS source matched for each object to avoid confusion with 
other sources in the field. For the purposes of probing the nuclear 
emission, we require that there are IRAS data available at least one 
of 12 or 25 /im (both are preferred, to get an impression of the 
shape of the IRAGN SED). 

We then enforce a quality criterion on the IRAS photometry, 
requiring at least one of the 12-100 /im points to have a signal-to- 
noise ratio (SNR) greater than 4.0, and not using any data with a 
SNR < 3.0. The photometry satisfying the above stringent criteria 
were used in further analyses. This yielded 71 out of 113 objects 
in the low-redshift sample. We identify the sources f or which pro- 
cesse d 9-month BAT data are available online fsee lTueller et al.l 
2008). We also require that an estimate of the black hole mass from 
the 2MASS Point Source Catalogu e (PSC) magnitudes is avail- 
able, using the method presented in IVasudevan et"aD (l2009h . The 
2MASS magnitude for the cD galaxy Cyg A yields a very low esti- 
mate of the black hole mass (log Mbh /Mq=1.9), since the bulge is 
unusually large in this particular galaxy . We therefore use the m ore 
reliable dynamical mass estimate from lTadhunter et aP (l2003h for 
this object (log Mbh/M0=9.4). This finally yields 64 objects for 
further analysis, listed in Table[T](29 low absorption and 35 high ab- 
sorption objects), constituting 57 per cent of the available low red- 
shift sample. Luminosity distances are computed from redshifts for 
most of the sample. Due to the large number of NGC galaxies in the 
sample, we use the redshift-independent distanc e estimates from 
the Nearby Galaxies Catalogue where available (iTullvll 19941) and 
scale luminosities appropriately. This only produces a pronounced 
change for two objects, NGC 4388 (where the change in luminosity 
from the redshift-estimated distance is a factor of ^ 4) and NGC 
4051 (where the luminosity changes by a factor of ^ 3). For all 
other NGC galaxies, the changes are negligible and less than the 
random errors from other sources. To allow easy comparison with 
previous work in this area and on the 9-month catalogue, we fi- 
nally bifurcate the sample on A^h, producing two subsamples of 
low (log Nh <22, 51 AGN) and high (log A^h >22, 61 AGN) 
absorption (using the dividin g line between 'ob scured' and 'unob- 
scured' sources suggested bv IWinter et al ]l2009h . 



3 GENERATING SEDS FROM BAT AND IRAS DATA 
3.1 BAT data 

The 4-channel BAT data are available readily processed online (see 
321) for analysis with the X-ray analysis software XSPEC. The key 
goal of using the BAT data is to obtain an accurate estimate of the 
intrinsic X-ray continuum. The BAT data become particularly use- 
ful at determining the intrinsic spectral shape and normalisation 
in the more obscured sources (provided they are not too Compton 
thick, Nh > lO^^cm-2). 

For our purposes, we wish to extrapolate the 2-10 keV lumi- 
nosities of these sources from their BAT data in order to calculate 
bolometric corrections. These luminosities also allo w comparison 
with th e Li2.3Aim — i>2-iokeV correlation presented in lGandhi et al.l 
(l2009h . The key assumption needed when extrapolating to calcu- 
late L2-iokeV is the spectral shape across the X-ray continuum, 
parameterised by the photon index F, assuming an intrinsic spec- 
trum consisting of a power-law of the form N(E) (x across 
the whole 0.1-200 keV range. The analysis of lTuelleretaLl ( l2008h 
suggests an average photon index of 2 in the BAT energy band 

^ http://swift.gsfc.nasa.gov/docs/swift/results/bs9mon/ 
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AGN 



A^H/(1022cm-2) 



redshift 



comments 



3C 382 

ESQ 548-G081 

Mrk 590 

Mrk 766 

Mrk 841 

NGC 4051 

NGC 7469 

NGC 985 

IRAS 05589+2828 

Mrk 79 

Mrk 279 

Mrk 509 

NGC 7213 

NGC 4593 

2MASX J21 140128+8204483 

NGC 5548 

ESO511-G030 

3C 390.3 

Mrk 290 

3C 120 

MCG -06-30-015 

MCG +08-11-011 

ESQ 490-G026 

NGC 3516 

NGC 931 

NGC 3783 

IC 4329A 

NGC 4945 

IRAS 09149-6206 

NGC 7314 

NGC 2992 

Mrk 3 

4U 1344-60 

NGC 526A 

NGC 3227 

EXO 055620-3820.2 

NGC 5506 

NGC 2110 

Mrk 6 

2MASX J04440903+28 13003 

NGC 6860 

ESQ 005-G004 

NGC 7582 

NGC 7172 

CygA 

PGC 13946 

Mrk 348 

Mrk 1498 

Mrk 18 

NGC 6300 

ESQ 103-035 

IC 5063 

NGC 4507 

NGC 4388 

MCG -03-34-064 

ESQ 297-018 



0.00+004 

0.0063 (upper limit) 

0.013 (upper limit) 

nip:+0.008 

U.U±0_Q QQg 

U.UZO_Q 

r. ^oi+O.Oll 

U.UC)i_Q 

047+0-023 

U.U^^ _Q Q21 

0.07tro5 

0.098t0:0^2l 
12+0-03 

n 15+0-03 

U.iO_o.05 

16+0-01 
19+003 

<;>P^ + 0.016 

^•^^-0.015 
n OO+0.04 
U.oo_o.o2 

0.353t0-32 
n S6+0-08 

n 57+0-21 
•^•^'-0.14 

61+0 03 
0.793t0-957, 

0.8510;?? 

1 16+0-01 

1 19+2-21 

1 45+0-20 
1 50+0-14 

1 74+0.12 

'^-0.09 
^ ^7+0.14 
-0.14 

2 84+0-19 

M.33 
-1.19 
o oQ+0.31 
<J-^i^_0.25 
M.33 
-1.30 

5 58+0-16 

^•^^-0.16 
f 1.46 
-1.00 

8 19+3-42 

°-^^-3.30 
11 0+21-0 

^±1 



3.26^ 



4.53^ 



7.39j 



16^ 
21 5+0-8 



21. 



612-6 

21.78+1^4 
34.28lt^° 
36.17l^J^ 

40 73+4-79 

41 71+4.70 



0.05787 radio galaxy 

0.01448 

0.026385 

0.012929 

0.036422 

0.002336 

0.016317 

0.043143 

0.033 

0.022189 

0.030451 

0.034397 

0.005839 

0.009 

0.084 radio galaxy 

0.017175 

0.022389 

0.0561 radio galaxy 
0.029577 

0.03301 radio galaxy 

0.007749 

0.020484 

0.02485 

0.008836 

0.016652 

0.00973 

0.016054 

0.001878 

0.0573 

0.004763 

0.00771 

0.013509 

0.012879 

0.019097 

0.003859 

0.03387 

0.006181 

0.007789 

0.018813 

0.011268 

0.014884 

0.006228 

0.005254 

0.008683 

0.056075 cD galaxy, radio galaxy 

0.036492 

0.015034 

0.0547 

0.011088 

0.003699 

0.013286 

0.011348 

0.011801 

0.008419 

0.016541 

0.0252 



continued on next page... 

Table 1. Objects used in the samp le. Values of absorbi ng column density A^h are taken from I Winter et al.l ( l2009h . * For NGC 4945, observations spanning 



2-20 keV from the Ginga satellite dlwasawa et al.ll 19931) yield a much higher value for A^h of 1024-'^cm 2^ placing it in the heavily 'obscured' category; the 
optical images reveal a heavily dust-extincted edge-on galaxy. 
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AGN 


A^H/(10^^cm~^) 


rcdshift 


connncnts 


3C 403 


45.01^;° 


0.059 


radio galaxy 


NGC 788 
ESQ 506-G027 


Ar on+4.68 

46.8914 47 
76.821^:^,^ 


0.013603 
0.025024 




NGC 1142 
NGC 5728 


79.75l^:gJ 
82.0i^;g 


0.028847 
0.009353 




NGC 3281 


+ 16.32 


0.010674 




NGC 1365 


104.8 


0.005457 




NGC 612 


129.701^^3^,° 


0.029771 


radio galaxy 



Table[T](continued) 



with a root-mean- square spread of ^ 0.3. We first consider the 
approach of calculating the intrinsic 2-10 keV luminosities by em- 
ploying r = 2 for all sources, but this yields systematically higher 
values of L2-iokey than those from X-ray analyses presented in 
IWinter et"aD (l2009h . Although the BAT 9-month data are averaged 
over many months and the X-ray analyses in Winter et al. ( 20ol) 
are from individual shorter observations, statistically one expects 
better agreement, assuming the X-ray analyses have managed to 
recover the true intrinsic luminosity. 

A second approach invo lves using the value s of the 2-10 keV 
photon indices presented in IWinter et al ] (l2009l) to constrain the 
shape of the overall X-ray spectrum, but using the BAT data to 
constrain the normalisation and hence the luminosity. The X-ray 
variability seen in the well-studied Seyfert 1 galaxy MCG-06-30- 
15 supports such an approach, as the observations are consistent 
with variation in no rmalisation while niaintaining a relatively con - 
stant spectral shape dVaughan & Fabianl2004ljMiniutti et al.l2007h . 
This produces better agreement and less of a systematic shift above 
the Winter et al. ( 2009) values, but still results in a large scatter. 
Part of this could be attributable to extreme values of F found for 
some of these sources in the X-ray analyses. On physical grounds, 
one expects a minimal photon index of 1.5 on canonical inverse- 
Comptonization scattering models for the coronal emission, and 
similarly, photon indices above 2.2 are likely to indicate com- 
plex absorption is preventing the recovery of the intrinsic shape, 
despite the inclusion of some absorption components. In our fi- 
nal approach, we therefore constrain V to He within the hard lim- 
its 1.5 < T < 2.2 when the values from I Winter et al. I (l2009l) lie 
outside this range. This correction is seen to be necessary for a 
number of the high absorpti on, 'complex' spectrum sources iden- 
tified in I Winter et al.l (l2009l) . The scatter between the BAT values 
and the Winter et al. values is significantly reduced, and the results 
of this comparison are presented in Fig. [T] We note that there are 
two objects which deviate significantly from the one-to-one cor- 
respondence line: these are NGC 4945 at logCL l^inlt'l)- 40.1 
(classified as unobscured in IWinter et"al ] l2009b and 3C 403 at 
log(L2^jQ^^Y)^ 42.8 (classified as obscured). Despite th e identifi- 
cation of the former as unobscured in lWinter et al ]l|200i), these are 
both known to have very complex spectra in which the 2-10 keV 
luminosity is prone to being underestimate d with simple absorbed 
power-law fits (see e.g. Iwas awa et al .Il993h . The BAT-extrapolated 
values for L2-iokeV are likely to be more representative in these 
cases. 

Having thus determined a sensible X-ray spectrum, we deter- 
mine the total X-ray luminosity by int egrating between .5 and 100 
keV. A cut-off of 500 keV (as used bv lPozzi et aDl2007h would not 
produce an appreciable difference in total X-ray luminosities for 
objects with F 2. The effect would be ^ 20 per cent in the X- 




log(L2.iokev) from Winter et al. (2009) 

Figure 1. Comparison between 2-lOkeV X-ray luminosity extrapolated 
from fitting to the BAT data with the values presented in Iwinter et"aD 
^2009 ). The photon indices (F) from Winter et al.,. |20Q^) were used in fit- 
ting the BAT data, but were restricted to He within 1.5 < F < 2.2 on phys- 
ical grounds. Filled circles and empty circles represent objects with high 
and low absorption respectively (with objects above log(A^H)=22 classed 
as having high absorption). 

ray luminosity, but since the dominant component of Lboi is the 
accretion disc emission (reprocessed or otherwise), this choice of 
cut-off does not have a big impact. Additionally, a 500 keV cut 
would be outside the BAT bandpass. Recent work by Mushotzky 
et al. (in prep) suggests that the majority of high-energy cut-offs in 
BAT spectra for the 22-month catalogue AGN are located at 100 
keV, which would support the use of the upper integration limit 
energy used here. 

We note that two objects in our high absorption sample could 
be borderline Compton-thick (A^h > lO^^cm"^) in I Winter et al.l 
(2009), namely NGC 612 and NGC 1365. In the case of NGC 
1365, this source has been identified to have variable absorption, 
with the source switching between Compton-thick a nd Compton- 
thin s tates on timescales of a few tens of kiloseconds jRisaliti et alj 
l2009h . There are additionally some sources which are identified 
as Compton-thick in the literature due to different ap proaches of 
modelling the spectr um (e.g. Mrk 3 jAwaki et aDl2008l NGC 4945, 
llwasawa et al.l 1993b . The g alaxy NGC 4945 in p articular may have 
iVn as high as lO^^'^cm"^ (llwasawa et al.ll993h . making it unsuit- 
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able for the type of Lboi calculation used here, and we therefore 
exclude it from our final bolometric correction results. The other 
very high-absorption objects represent the limits of the approach 
outlined here and their results must be treated with caution. 



3.2 IRAS data 

The next step in the process is to determine the total nuclear IR 
luminosity, which, under the reprocessing paradigm, constitutes 
the remaining major component of the bolometric luminosity. The 
IRAS 12 and 25 /im fluxes are thought to be dominated by the nu- 
clear component, but they are typically calculated for apertures sig- 
nificantly larger than the expected torus dimensions - indeed for the 
vast majority of local AGN, the torus cannot be resolved even with 
very high resolution imaging. There is therefore a likelihood of sig- 
nificant non-nuclear flux contaminating the photometry, which we 
now turn to address. 



3.2. 1 Host galaxy/starburst contamination in IRAS fluxes ( I ): 
using the Li2/xm — i^2-iokeV relation 

iGandhi et all (l2009h present a strong correlation between nuclear 
12.3 fim and 2-10 keV luminosities for AGN. Assuming the colour 
correction between 12 /im and 12.3 /im is negligible, we can com- 
pare the relationship between ouvIRAS 12 /x m luminosity and BAT - 
derived 2-10 keV luminosity with that from lGandhi et al.l (l2009h , 
to get a statistcal estimate of the degree of host galaxy contamina- 
tion present in the IRAS fluxes at different X-ray luminosities. We 
present these data in Fig.|2l in a ddition to the co rrelation found for 
well-resolved sources from iGandhi et al.l (l2009h . This comparison 
also assumes that the flux discrepancies introduced by the different 
filters on IRAS and the VLT are negligible (the broad IRAS 12 fim 
filter could include many features contaminating the continuum, 
whereas the much narrower 12.3 fim VLT filter is likely to provide 
a better sample of the continuum). 

It is clear that the non-nuclear excess in the 12 /xm luminos- 
ity increases with lower X-ray luminosity; this is expected since at 
lower intrinsic nuclear luminosities (as traced by iv2-iokev), it is 
more difficult for the nucleus to outshine the host galaxy. This ex- 
cess is likely to be due to star-formation . We determine the excess 
with respect to the iGandhi et all (l2009h relation for all the AGN 
with detections in IRAS and plot them against L2-iokeV in Fig. [3] 
We also include upper limits for 39 AGN without IRAS detections, 
taking care to check whether the objects in this list were in a part 
of the sky surveyed by IRAS. We take the completeness limit of the 
IRAS FSC to be the upper limiting flux for these objects, ^ 0.2 Jy 

B 

We use the Astronomy Survival Analysis (ASURV) package 
from the StatCodes suite of utilities developed by Eric Feigelson 
Q to determine the correlation between the excess Alog(Li2yum) 
and logL2-iokeV including the effects of the upper limits. The 
'EM' and 'Buckley- James' algorithms available in ASURV yield al- 
most identical results. We obtain the following non-nuclear 12 fim. 
excess: 

Alog(Li2^m) = 19.988 - 0.4578 log(L2-iokev). (1) 




log(L2.iokev) from fit to BAT data 

Figure 2. Infrared 12 /im luminosity against 2-lOkeV luminosity from 
fits to BAT data. The solid line represent s the correlation found for well- 
resolved sources from lGandhi et all (l2009l) . Filled circles and empty circles 
represent objects with high and low absorption respectively (with objects 
above log(A^H)=22 classed as having high absorption). 



1.51 

♦ 




41 42 43 44 45 

log(L2.iokev) 

Figure 3. Offset AlogLi2^m between IRAS 12 /im fluxes and the relation 
of Gandhi et al. ( 2009), plotted against i^2-iokeV- Detected objects are 
represented using black filled circles (obscured objects) and empty circles 
(unobscured objects), and upper limits are shown using downward pointing 
arrows. The solid line is the best- fit taking upper limits into account (using 
the EM algorithm in the ASURV utility), and the dashed line shows the fit 
from ignoring non-detections. Over the luminosity range of the sample, the 
differences are not significant. 



^ http://irsa.ipac.caltech.edu/IRASdocs/surveys/fsc.html 
^ http://astrostatistics.psu.edu/statcodes/sc_censor.html 
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Figure 4. Infrared 12 /im luminosity against 2-lOkeV luminosity from fits 
to BAT data, corrected for the excess starburst/host galaxy flux seen in the 
large apertures for IRAS. The excess host contribut ion was obtained b y 
comparing the correlation from our data with that from lGandhi et alj l2009h . 



The correction exhibits a significant dependence on luminos- 
ity: for AGN with log(L2-iokev) ~ 41, the total IR 12 /xm flux is 
estimated to be 20 times larger than the nuclear flux alone, re- 
ducing to zero at X-ray luminosities of log(L2-iokev) ^ 43.5. We 
then correct the 12 /im luminosities according to this correction (up 
until the threshold at which the correction becomes negative) and 
re-plot the results in Fig.|4] 

In this approach, we assume that the 25 /xm luminosity needs 
to be corrected by the same factor, despite a possible larger non- 
nuclear contribution at longer wavelengths. The larger excess ex- 
pected at 25 /im would further decrease the total IR nuclear lumi- 
nosities obtained. However, since we do not have any information 
on the 25/im-L2-iokeV correlation, an estimation of the extra cor- 
rection required is beyond the scope of this paper. 

The Li2/xm — i^2-iokeV corrclation in Seyferts itself presents 
a possibility for predicting the approximate range of bolometric 
corrections under the reprocessing paradigm, since the Li2/xm — 
^2-iokeV relation links the two bands used in our calculation of 

^bol- 

We take the bolometric luminosity Lboi to be the sum of the 
total integrated IR luminosity and the total X-ray luminosity, i.e. 
Lboi = ^iR,tot + ^x,tot. Assuming that the monochromatic 12 
/im luminosity scales up to the total integrated luminosity by some 
factor K, dependent on the IR SED and including the effects of 
geometry and anisotropy (see ^3.2.31 for details), we can re- write 
this as: 



K Ll2/xm + i>x,t- 



(2) 



This can be converted into an expression for the bolometric correc- 
tion A>:2-iokeV by dividing through by L2-iokev: 



The final part of the expression Lx,tot/i>2-iokeV is constrained to 
lie within 3.1 and 5.3 by the range of photon indices 1.5 < F < 2.2 
adopted here as discussed in ^3.11 We can then include the Li2/xm — 
^2-iokeV relation as Li2^m — 10'^iv2_iokev to finally obtain 



At2-10keV = K 10"iv2^-LQ^j^gy + -— 

^2 



(4) 



^2-lOkeV 



K- 



^12/xm 



^X,tot 



(3) 



The distribution of bolometric corrections seen will therefore de- 
pend on the diversity of IR SED shapes and assumptions about 
torus geometry and emission (via K), the degree of spread in- 
trinsic to the measured Li2/xm — i^2-iokeV relation (via a and 
h) and the amount of variation in the X-ray spectral shape (via 
^x,tot /^2-iokev)- The nuclear SED templates from ISilva et al.! 
(I2OO4I) used in this paper, in conjunction with assumptions about 
torus geometry and anisotropy of emission ( ^3.2.31 l yield a range in 
of 4 — 15. Comb ining this with the spread in a and h reported in 
iGandhi et al.l (l2009h and the range of X-ray SED shapes discussed 
above, this predicts bolometric corrections ranging from 6-60 over 
the luminosity range 41.0 < log(L2-iokev) < 45.0 spanned by 
the objects considered here, with an increase in Ai:2-iokeV with 
^2-iokeV as predicted by equation Therefore, by requiring that 
our objects lie on the Li2/xm — i>2-iokeV relation, we naturally ex- 
pect a certain distribution of bolometric corrections, with the de- 
gree of spread related to the tightness of the Li2/xm — i>2-iokeV 
relation observed. The tighter the distribution, the smaller the range 
of bolometric corrections will be, centering around 10-30. Since, 
under the assumptions used in our paper, the requirement of a 
tight Li2/xm — i^2-iokeV relation can limit the range of bolo- 
metric corrections allowed, we also consider an approach for re- 
moval of the non-nuclear component which is independent of the 
i^i2/xm — iv2-iokeV relation. 



3. 2. 2 Host galaxy/ starhur St contamination in IRAS fluxes (2): 
using host galaxy IR SED templates 

We also employ the strategy used bv lPozzi et al.! (l2007h to account 
for the host gal axy, namely fitting the nuclear and host galaxy SED 
templates from lSilva et alJ (l2004h to the broad-band (12-100 /im) 
IR data. This does not presuppose any relation between Lir and 
^2-iokeV, but we can also determine whether this approach to cor- 
recting for the host galaxy reproduces the Li2^m — i^2-iokeV re- 
lation. We select from the h ost galaxy S ED templates for different 
luminosity bins presented in lSilva et al.l ([2004) using our estimate 
of L2-iokeV from the BAT data. In a few cases where 60 or 100 
/im data were not available, the normalisation of the host template 
was tied to the nuclear template normalisation, to provide a sen- 
sible estimate of host galaxy contamination based on the observa- 
tions of Silva et al. (2004) in their sample of AGN. Some example 
SEDs with both host and nuclear SED templates fitted are shown in 
Fig.[8] We also present the distribution of Li2/xm against L2-iokeV 
in Fig. [21 with the nuclear Li2/xm determined from the 12 /im flux 
of the nuclear part of the SED fit only. 

In the absence of high-resolution data for a large fraction of 
our sample, these methods for removing the host galaxy can only 
go so far in tackling this complex issue. In Fig.[6]we present a com- 
parison of our uncorrec ted IRAS fluxes wit h the accurate nuclear 
VLT/VISIR fluxes from lGandhi et al.l (l2009l) . It can be seen that in 
some cases, the degree of host contamination in the IRAS fluxes can 
be very large (e.g. NGC 4388); while in others the IRAS fluxes are 
almost identical to the nuclear fluxes from VISIR/VLT (e.g. NGC 
6300), indicating that no correction for contamination is necessary. 
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3.2.3 Total nuclear IR luminosity 



42 43 44 

log(L2.iokev) from fit to BAT data 

Figure 5. Nuclear 12 /xm luminosity against 2-lOkeV luminosity from fits 
to BAT data, from combined host and nuclear IR SED template fitting to 
the IRAS 12-100 fim data. 
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GC 7314 

eNGC 4388 



41. 5L 



42 42.5 43 43.5 44 

log(L2_io) from fit to BAT data 



44.5 



Figure 6. Infrared 12 /im luminosity against 2-lOkeV luminosity from fits 
to BAT data. Black points are IRAS data uncorrected for host galaxy con- 
tamination, and empty circles show the nuclear 12 yum luminosity from the 
VISIIWLT data presented in .Gandhi et al.. C2009.) . 



; way m ac 

aperture effects, but as discussed bv lHorst et al.l (|2008i), higher res- 
olution IR data which isolates the nuclear emission is ultimately 
much preferred. 



We broadly follow the approach of IPozzi et al ] (l2007h outlined 
in their study on high-redshift, luminous, obs cured quasars and 



fit the nuclear and host galaxy templates from ISilva et ah I (l2004l) 



to our data. The appropriate nuclear template from Silva et all 
(l2004l) is selected b ased on the absorbing column A^h reported in 
IWinter et"aD (l2009h . The fitting is performed using a simple least- 
squares algorithm from the SCIPY suite of functions available for 
use with the PYTHON programming language. If the host galaxy 
contamination is removed via method 1 ( ^3.2.11 ) above, we only fit 
the 12-25 /xm data and correct the integrated nuclear IR luminosity 
by the fraction given by equation If method 2 ( ^3.2.21 ) is used, 
only the nuclear component of the fit is integrated. The observed IR 
luminosity is then obtained by integration of the nuclear template 
between 1 and 1000 /xm. 

We reiterate briefly here that the absorptions seen in X-rays 
and the amouts of dust responsible for reddening the optical-UV 
SED/produc ing the reprocessed IR do not necessarily match (e.g. 
lMaiolino"& Risalitl 2007i) . This calls into ques tion the use of th e 
X-ray A^h for selecti ng the template to use from lSilva et al.l(l2004l) : 
however, the fact that lSilva et al.l (|2004|) find from observations that 
the IR SED shapes of AGN can be broadly grouped on A^h would 
argue for at least a statistical correlation between the two types of 
absorption, in addition to the findings from surveys discussed in ^ 
In this first-order approach, the selection of nuclear SED template 
based on A^h is therefore sensible. 

Under method (2) for host galaxy removal, a few objects (Mrk 
841, NGC 612, PGC 13946, 2MASX J04440903+28 13003, EXO 
055620-3820.2, 4U 1344-60, ESQ 297-018, NGC 7314) yield a fit 
with zero nuclear component and only a host galaxy component. 
These probably represent cases where the host galaxy contamina- 
tion is very large and the nucleus is too buried to get an accurate 
estimate of its presence from our simple SED fitting alone. We ex- 
clude these objects from the results from method (2) due to these 
uncertainties. A detailed exploration of the different host and nu- 
clear SED models may be able to yield more satisfatory fits in these 
cases, but is beyond the scope of this paper. 

To convert this observed IR luminosity into a measure of 
the reprocessed nuclear accretion disc luminosity, correction fac- 
tors are required to account for the geometry of the torus and the 
anisotropy of the torus emission. We adopt the same approach as 
IPozzi et ^ (l2007h in correcting for these effects. The correction 
for geometry relates to the covering factor / which obscures the 
optical-UV emission from the accretion disc from view. Their ge- 
ometry correction is based on a statistical argument, employing the 
ratio of obscured to unobscure d quasars as foun d by recent X-ray 
background synthesis models (iGilli et al.ll2007b to infer a typical 
torus covering fraction of / ^ 0.67. This value is also consistent 
with the cove ring fraction obtained from recent detailed clumpy 
torus models jNenkova et aLll200"8bl) , for typical torus parameters 
(e.g. number of line-of-sight clouds ~ 5, with an opening angle of 
^ 30 — 45°). Inverting this, we obtain a factor ^ 1.5 by which 
the observed IR nuclear luminosity needs t o be multipl i ed to o btain 
the total optical-UV accretion emission. IPozzi et al ] (l2007h esti- 
mate the anisotropy correction factors by computing the ratio of 
the luminosities from face-on vs. edge-on AGN (with the column 
density A^ h param eterising the inclination of the torus). They nor- 
malise the lSilva e t al. (200^ templates for all different A^h values 
to have the same 30-100 /xm luminosities. The luminosity ratios 
for face-on to edge-on AGN are then calculated in the 1-30 /xm 
range, where the effects of anisotropy are most pronounced. They 
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obtain values of - 1.2 - 1.3 for 10^^ < iVn < lO^^cm-^ and 
^ 3 — 4 for A^H > lO^^cm"^ sources. For face on (low absorption 
logNn < 22) sources, no correction is necessary. We adopt a sim- 
ple approach, using values of 1.3 for 22 < logNn < 24 sources 
and 3.5 for logNn > 24 sources. 



4 CALCULATING THE TOTAL POWER OUTPUT 

We present some of the SEDs constructed from the IRAS and BAT 
data in Figs. |7] and [8] The bolometric luminosity Lboi is calcu- 
lated as the sum of the nuclear IR luminosity LiR,corr (corrected for 
torus geometry, anisotropy of emission and non-nuclear contamina- 
tion as detailed above) and the total X-ray luminosity, Lo.5-iookeV • 
For the high-absorption objects, we multiply the power-law by an 
intrinsic absorption comp onent (the WABS m odel in XSPEC, us- 
ing the value of A^h from I Winter et al.l (l2009h ) before calculating 
^0.5-iookeV to avoid double-counting part of the X-ray luminos- 
ity which i s reprocessed to the IR, in line with the approach of 
IPozzi et al] ([2002). For low-absorption objects, the difference be- 
tween the intrinsic and absorption-corrected 0.5-100 keV lumi- 
nosities are found to be negligible, so we leave out this step for 
that class. We also extract the absorption-corrected L2-iokeV for 
calculating bolometric corrections /^2-iokeV = i>boi/^2-iokeV. 

We also calculate Eddington ratios for our sample using black 
hole mass (Mbh) estimates. For all but one object (Cyg A), 
these are calculated from the Mbh — ^k, bulge relation (where 
^buige is the host galaxy K-band bulge l u minos ity), for con- 
sistency with the work of IVasudevan et al ] (I2OO9I) . We use the 
Im arconi &Huntl(l2003h formulation for this relation, using an iden- 
tical method to determ ine the bulge luminosity as that outlined in 
IVas udevan et al.' ('2009*). For Cyg A, the dynamical mass estimate 
from Tadhu nter et al.. C2003i) is used as discussed in 321 Eddington 
ratios are calculated using AEdd = ^boi/i^Edd for Eddington lu- 
minosities LEdd = 1.3 X 1O^^(MBH/M0)ergs"\ 

4.1 Calibrating the method 

We first attempt to estimate the accuracy of our bolomet- 
ric luminosity calculat i on using IR and har d X-ray emission. 
iBrocksoDP etaP (l2006h . IVasudevan & FabianI ([2009») (and refer- 
ences therein) outline the approach of using simultaneous optical, 
UV and X-ray data for calculation of intrinsic SEDs, and this ap- 
proach is developed fu rther using Swift XRT and UVOT data in 
IVasudevan et al.l J2OO9I) for low absorption objects in the Swift/BAT 
9-month catalogue. Using IR along with long-exposure BAT data 
produces a long-term averaged estimate of Lboi in contrast to 
the 'snapshot' approach from using simultaneous optical-to-X-ray 
data, but statistically one expects reasonable agreement between 
the two approaches, since the effects of variability should be av- 
eraged out when considering a large enough sample. We present 
comparisons be tween the value s of L boi, Ai:2-iokeV, AEdd and 
L2-iokeV from Vasud evan et al ] (l2009h and this study in Figs. [9] 
and[TOl for the two different host-contamination removal methods. 
There are 17 objects overlapping between the two studies. 

The bolometric luminosities determined from both methods 
agree reasonably well with those determined from UVOT and XRT 
fits, but both methods show a systematically larger Lboi than that 
determined from UVOT and XRT data. In particular, the degree of 
host galaxy contamination removed from fitting nuclear and host 
SEDs (method 2) yields estimates for Lboi about 0.3-0.4 dex larger 
than the optical-to-X-ray estimates. This could be attributed to a 



number of factors; firstly it is possible that the simple SED fitting 
(in method 2) underestimates the host galaxy contribution and more 
correction is necessary; secondly the geometry and anisotropy cor- 
rections assumed in scaling the observed IR luminosity could be 
too large for many of the objects where poor agreement is re- 
ported, and thirdly, the estimates from UVOT-XRT could them- 
selves be too small (the latter two phenomena would apply equally 
to method 1 and method 2). The Lboi from the UVOT and XRT fit 
depends on the black hole mass and tends to increase as the black 
hole mass decreases: IVasudevan et"aD (l2009l) discuss the possibil- 
ity that the mass could be overestimated by about a factor of ^ 2 
which could alleviate the discrepancy seen here. The problem is a 
complex one and ultimately requires better knowledge of the torus 
covering fraction (which ma y have a very broad distribution, e.g. 
iRowan-Robinson et al.ll2009l) and emission characteristics to per- 
form this comparison more accurately. 

The distributions of bolometric correction and Eddington ra- 
tio are more clustered, so a correlation between the two approaches 
is difficult t o discern. The finding o f a low distribution of Edding- 
ton ratios in lVasudevan et al.l (l2009h is, however, broadly confirmed 
with the IRAS and BAT method. The notable outlier Mrk 590 is 
known from the literature to have a very reddened optical-UV 
spectrum (confirmed with UVOT), and constraining the big blue 
bump shape using the K-band mass estimate yields a poor fit in 
I Vasudevan et al.l i2009h . If for this particular source, we employ 
the results obtained from using the reverberation mass instead, its 
accretion rate lies much closer to the line of one-to-one corespon- 
dence, moving to join the cluster of points around AEdd 0.02. 
Despite the scatter however, the general trend for agreement be- 
tween the two methods motivates us to continue in our attempts at 
addressing the key questions of this study. 

We also present some of the SEDs for the objects with XRT 
and UVOT data, now including the IRAS and BAT data (Fig.fTTT). 
There is some variation between the photon inde x used for fitting 
the BA T data and that found from the XRT data in I Vasudevan et al.l 
(l2009h . but generally both datasets show similar values. 



5 RESULTS AND DISCUSSION: BOLOMETRIC 
CORRECTIONS AND EDDINGTON RATIOS 

We present the bolometric corrections and Eddington ratios ob- 
tained in Figs. [12] and [13] The distributions are plotted as his- 
tograms in Fig. [14] and numerical values are provided in Table [3 
The values obtained from fitting a nuclear SED template to the IR 
data without any host galaxy correction are shown using empty cir- 
cles for comparison, and objects from the 3C catalogues (3C 382, 
3C 120, 3C 390.3 and 3C 403) are higWighted with blue squares, 
for easy identification in case some degree of jet contribution to the 
flux affects the results for these objects. 

For low-absorption objects, the statistical correction method 
for removing the host galaxy results in bolometric corrections of 
19 on average with Eddington ratios centered around ^ 0.035 
(Fig. [T2I left panel). In contrast, using the SED fitting approach 
yields (A>:2-iokev) ^ 24 with (AEdd) ^ 0.039 (Fig.[T3l left panel). 
Both are generally consistent with the fi nding of low bolometri c 
corrections for unabsorbed objects from I Vasudevan et al.l (l2009h . 
but the host SED fitting approach does have a greater proportion of 
objects at higher bolometric corrections. 

The object consistently exhibiting the largest bolometric cor- 
rection in the high absorption subsample is NGC 1365, which is 
one of two possible borderline Compton-thick objects in that sam- 
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Figure 7. Some example SEDs with IRAS and BAT data. Data points used in the fit (IRAS and BAT) are shown with filled black circles. Model fits are 
represented by the solid black lines. Only 12-25 /im IRAS data are used to fit the IR nuclear templates; longer wavelength (lower energy) 60 and 100 /im data 
are shown for information (unfilled circles). The nuclear luminosities from the fit are corrected for host contribution using method 1 ( ^3.2. U . 




Figure 8. Some example SEDs with IRAS and BAT data, with both host and nuclear SED fitting ( ^3.2.2t . Data points (IRAS and BAT) are shown with filled 
black circles. The full model fit is represented by solid black lines; the dotted and dot-dashed lines show the host and nuclear components respectively. 



pie and known to have variable absorption as discussed in ^3.11 If 
the source was predominantly in a Compton-thick state during the 
BAT period of observation, the BAT flux would be substantially 
lower than the intrinsic AGN flux as pointed out in ^3.11 giving an 
artifically high bolometric correction. Additionally, this object is 



known to exhibit a nuclear starburst JStrateva & Komossall2009l and 
references therein) so is particularly susceptible to a high level of 
starburst contamination in the IR fluxes, which may not have been 
completely removed by either method of host galaxy removal. Both 
of these effects would serve to increase the bolometric correction to 



© 0000 RAS, MNRAS 000, 000-000 



The power output of local obscured and unobscured AGN 1 1 




?.Edd (UVOT+XRT) log(L2.iokev/erg s"^) (XRT) 

Figure 9. Comparison between values of bolometric luminosity, 2-lOkeV bolometric corrections, Eddington ratios and 2-lOkeV luminosities using the time- 
averaged, non-contemporaneous IRAS and BAT data, and the method from lVasudevan et all (l2009h using simultaneous optical, UV and X-ray data from Swift. 
The solid lines represent the desired 1:1 correlation, and empty circles show the results before correcting for host galaxy contamination of IRAS fluxes. Host 
galaxy correction to //M5-determined IR luminosities as detailed in ^3.2.11 



an abnormally high value. A more detailed object- specific analysis 
is needed to check whether the bolometric correction for this object 
is intrinsically high. 

For high-absorption objects, the statistical approach to cor- 
recting for the host galaxy (Fig. [TJl right panel) yields a wider 
spread of bolometric corrections than for low-absorption objects, 
but values are still centred on low values ((Ai:2-iokev) ^ 15) with a 
predominantly low Eddington ratio distribution ((AEdd) 0.034). 
However, under the host SED correction method, some notable 
outliers with negligible host galaxy component are located near 
AEdd 0.1 with bolometric corrections of ^ 80. These outliers 
are (from highest bolometric correction downwards) NGC 1365, 
IC 5063, MCG -03-34-064 and NGC 3281. In the case of NGC 
3281, the degree of host contamination is known to be small from 
the comparison with VLT/VISIR fluxes in Fig.[6l so the high accre- 
tion rate and value of K2-iokeV may be real for such objects (bear- 
ing in mind the caveats regarding NGC 1365 discussed above). In- 
spection of the host and nucleus SED fits for these four objects in 
Fig. [15] implies sensible estimation of the host contamination, im- 



plying that method 1 has overestimated the nuclear contamination 
for these sources. The large bolometric corrections from method 2 
for these objects may therefore, on average, be closer to the true val- 
ues. These objects could then potentially be part of the 'transition 
region' to hig h bolometric corrections at and above AEdd ^ 0.1 
postulated bv lVasudevan & FabianI (120071) . (2009). Under method 
2 for host galaxy correction, at AEdd < 0.03, bolometric cor- 
rections have an average value of ^ 15, increasing to a value of 
^ 32 for 0.03 < AEdd < 0.2 (the average Eddington ratio for 
the whole sample is 0.043). This could be a preliminary in- 
dication of a significant minority of higher- accretion rate objects 
amongst high-absorption AGN and may provide tentative evidence 
th at the Eddington ratio-depe ndent bolometric correction scheme 
o^Vasude van & FabianI ( [2007h may be borne out in high-absorption 
AGN. Better constraints on the nuclear MIR emission and black 
hole mass are needed to confirm this, however. 

Under the statistical method for host-galaxy correction 
(method 1), fitting both 25 and 12 /im points for objects with par- 
ticularly 'red' 12-25 /xm SEDs to a nuclear SED template results 
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Figure 10. Comparison between values for bolometric luminosity, 2-l OkeV bolometric correc tions, Eddington ratios and 2-lOkeV luminosities from IRAS 
and BAT data with the results from Swift UVOT-XRT data presented in I Vasudevan et al.1 i2009h . Host galaxy in IR accounted for as detailed in ^3.2.21 Key as 
in Fig. [21 



in a fit inconsistent with the 12 /xm datum. Such SED shapes prob- 
ably imply that the nuclear 25/im luminosity is primarily from the 
host-galaxy rather than the nucleus, and fitting to only the 12/im 
point would be more robust. We estimated the effect of this on the 
results by fitting only the 12 /im for the handful objects with such 
SEDs, and find that the resultant bolometric luminosities are altered 
on average by a few per cent (the largest change is for NGC 7582, 
causing a reduction of 40 per cent in the total IR luminosity). Over- 
all, such modifications do not produce any significant changes in 
the results. More detailed model fits to better quality IR data in the 
future could be fit using a more considered approach. 

We note that a few AGN excluded from the 'host plus nucleus' 
SED fit subsample due to the fit yielding zero nuclear contribution 
(see ^3.2.31 ) are more likely to have very small nuclear IR luminosi- 
ties, therefore would be expected to lie at lower bolometric correc- 
tions and Eddington ratios; this would not be expected to alter the 
distributions seen significantly. 

These results are illuminating for our understanding of accre- 
tion in local AGN. Firstly, they confirm that the majority of the 
Swift/BAT catalogue AGN are accreting at low Eddington ratios 



dWinter et al.ll2009l) across the range of available absorption prop- 
erties probed. This result is not significantly modified by possible 
over-estimates of the black hole n iass using our Mbh ~ ^ buige ap- 
proach, as discussed in detail in I Vasudevan et al.l (l2009l) , since if 
we calibrate this approach against reverberation mapping as done 
in their paper, we obtain an offset of a factor of 2, which would 
shift the Eddington ratios obtained here to correspondingly higher 
values. However, the centre of the distribution would still be lo- 
cated at low values. At any rate, there are geometrical uncertainties 
in reverberation mapping which give rise to a tolerance of a factor 
of 3 in calibrating the mass estimates. 

However, we note the presence of five AGN for which the Ed- 
dington ratios co i ild be uncertain because of their mass estimates. 
IVasudevan et"aD (l2009h use the 2MASS Point Source Catalogue 
(PSC) flux as an estimate of the total unresolved nuclear and bulge 
light, and calculate the fraction of the 2MASS PSC flux which can 
be attributed to the bulge. Inspection of the 2MASS images for 
all of the nearby galaxies reveal that there may be a few sources 
for which the assumption of an unresolved bulge is not valid, and 
(apart from the case of Cyg A already discussed) those overlapping 
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Figure 11. A selection of SEDs with IRAS, UVOT, XRT and BAT data. Key as in Fig.|7]with the following additions: crosses represent simultaneous UVOT 
and XRT data, heavy dotted lines represent the fits to those data with intrinsic absorption included (optical-UV redde ning and X-ray absorpti on) and dashed 
lines represent the de-reddened, absorption corrected SED fit from which bolometric luminosities were estimated in Vasudeva n et all l2009h . The top three 
SEDs show the fitting of only a nuclear IR SED with host galaxy contamination corrected for afterwards as detailed in ^3.2.11 the lower three show the 
combined fitting of nuclear and host galaxy SED templates in the IR as described in ^3.2.21 




100 



10 







' ' ' ' 1 

o 


log Nh > 22 


o 




— i 


o 

~k~ o o 


—i 








' o ; 
1 1 1 1 1 



10-^ 



0.01 



0.1 



Figure 12. Bolometric corrections against Eddington ratios for low absorption (left panel) and high absorption (right panel) objects, with statistical correction 
for non-nuclear contamination of IRAS fluxes (black filled points) as detailed in ^3.2.11 Results before applying the contamination correction are shown for 
comparison (empty circles). Objects highlighted with blue squares are from the 3C catalogue; these radio-loud objects are demarcated to identify objects in 
which some jet contribution might be affecting the results. 
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Table 2. Table of results from analysis of IRAS-BAT SEDs. Luminosities L are presented as l og(L) with L given in ergs~^. a - Total (0.5-100 keV, first 
value), followed by 2-10 keV (second value) X-ray luminosity, using the photon index F from I Winter et al.1 ( l2009h (requiring 1.5 < F < 2.2). Errors (in 
luminosity) are typically ~ 10 — 20 per cent, b - total integrated nuclear IR luminosity, with geometry, anisotropy and non-nuclear contamination correction 
factors applied (see text). The first value includes the statistical host galaxy correction based on i^2-iokeV' ^^d the second value is that obtained from combined 
nuclear and host SED template fitting. Errors are typically ~ 5 — 10 per cent, d - Bolometric luminosity (i^o.5-iookeV + ^i-iooo/xm)- Values from both 
methods of host galaxy removal provided, as for c. Errors are typically ~ 5 per cent, d - 2-lOkeV bolometric correction ^C2-iokeV = ^boi/^2-iokeV' 
with statistical host galaxy removal/host-i-nucleus SED fitting methods, e - Bolometric correction for observed nuclear 12/im monochromatic luminosity 
= Lboi/(AL^^^™), for two different methods of host galaxy removal. / - l ogarithm of black hole mass from K-band bulge magnitude estimate (in 
solar masses). Errors are typically ~ 0.1 dex, but for a discussion of systematics see I Vasudevan et al] l2009h . * Masses mar ked with an asterisk are likely to 
have the bulge luminosity, and hence black hole mass, underestimated: see text for details, f We use the dynamical mass from lTadhunter et all l2003h for Cyg 
A; since its bulge is abnormally large by virtue of being a cD galaxy, the use of the 2MASS PSC is not valid and would underestimate the mass by a factor of 
~ 30. - Eddington ratio Asdd = ^boi/^Edd^ using both methods of host galaxy removal. Random errors are typically ~ 10 — 20 per cent. 



© 0000 RAS, MNRAS 000, 000-000 



The power output of local obscured and unobscured AGN 15 




^Edd ^Edd 

Figure 13. Bolometric corrections against Eddington ratios for low absorption (left panel) and high absorption (right panel) objects, where the host 
galaxy/starburst contribution in the IR has been removed by SED fitting as detailed in ^3.2.21 key as in Fig. [El 



AGN 






Li 


-lOOO/xm^ 












Mbh-^ 


AEdd^ 


NGC 3281 


43.3,42 


8 


44 


1/44 


4 


44 


2/44 


4 


21.4/40.1 


6.38/4.57 


7.20 


0.071/0.13 


NGC 4388 


43.0,42 


6 


43 


6/43 


6 


43 


7/43 


7 


11.4/12.8 


7.83/5.26 


7.07 


0.032/0.035 


NGC 4507 


43.8,43 


3 


44 


1/44 


2 


44 


3/44 


3 


10.3/11.6 


8.17/5.83 


7.70 


0.030/0.034 


NGC 526A 


43.7,43 


3 


44 


2/44 


3 


44 


3/44 


4 


12.0/14.2 


6.79/5.87 


7.93 


0.020/0.023 


NGC 5506 


43.5,43 





44 


0/44 


3 


44 


1/44 


3 


12.3/21.3 


7.80/5.53 


7.67 


0.020/0.035 


NGC 5728 


43.3,42 


7 


43 


5/43 


1 


43 


7/43 


5 


9.70/6.40 


10.3/10.2 


7.15* 


0.028/0.018 


NGC 612 


43.7,43 


7 


44 


4/- 




44 


5/- 




6.24/- 


3.01/- 


8.01 


0.022/- 


NGC 6300 


42.4,42 


1 


42 


7/42 


6 


42 


9/42 


8 


7.06/5.67 


8.26/7.17 


6.70* 


0.013/0.0095 


NGC 6860 


43.3,42 


6 


43 


6/44 





43 


8/44 


1 


14.6/25.6 


6.25/5.63 


7.91 


0.0059/0.011 


NGC 7172 


43.4,42 


9 


43 


6/43 


5 


43 


8/43 


7 


8.16/7.66 


7.53/7.94 


7.36 


0.020/0.018 


NGC 7314 


42.7,42 


4 


42 


71- 




43 


0/- 




3.86/- 


11.3/- 


6.14 


0.053/- 


NGC 7582 


42.5,41 


8 


43 


3/43 


3 


43 


3/43 


3 


34.4/33.7 


9.23/5.53 


7.44* 


0.0060/0.0058 


NGC 788 


43.4,43 


1 


43 


7/43 


9 


43 


9/44 





5.50/8.06 


6.67/5.56 


7.39 


0.023/0.034 


PGC 13946 


43.9,43 


3 


44 


21- 




44 


4/- 




12.1/- 


-/- 


7.46 


0.064/- 



Table 12] (continued) 



with our sample are NGC 6300, NGC 3227, NGC 5728, NGC 7582 
and NGC 1365. If the bulge luminosity has been underestimated in 
these sources, the small resultant black hole masses will artificially 
increase their Eddington ratios. If corrected, this should shift the 
average Eddington ratio down for this subsample. 

We note that the bolometric luminosities are not dependent on 
the mass in the simple approach used he re, in contrast to the method 
outlined in lVasudevan & FabianI (l2007h and their subsequent stud- 
ies, which involves fitting a multicolour accretion disc model to 
optical-UV data with the model normalisation constrained by the 
black hole mass. In the zeroth-order approach presented here, the 
main uncertainties in determining Lboi are the accuracy of the cor- 
rection for non-nuclear contamination, the accuracy of the nuclear 
IR SED templates used and the degree of hard X-ray variability 
JPozzi et a l. 2007 find that using various different IR SED tem- 
plates does not significantly alter the extrapolated bolometric lumi- 
nosity, however). As a result, uncertainties in mass estimates only 
affect Eddington ratios, not bolometric luminosities. 



Bolometric corrections are gene rally lower than those ob- 
tained for high Ed dington-rate sourc es dVasudevan & Fabianll2009l . 
2007) and quasars (lElvis et all 19 94). These results provides a very 
interesting window onto the properties of the more complex high 
absorption AGN, as hinted at previously bv lPozzi etal1(l2007h . Our 
approach aims to by-pass the unwanted complexities introduced by 
absorption in the optical-to-X-ray regime by using the reprocessed 
IR and (relatively) unaffected hard X-ray BAT data. Indications 
of predominantly low bolometric corrections for high absorption 
AGN have consequences for matching the local SMBH mass den- 
sity with that inferred from the X-ray background. A significant 
population of highly a bsorbed AGN is required to fit the X-ra y 
background spectrum ('Gandhi & FabianI I2003L iGilli et all 120071) . 
and as discussed in iFabian (.2004b . bolometric corrections of around 
^ 10 — 20 are appropriate for reconciling the XRB with the local 
black hole mass density. 
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2.51 I 2.5 




log(XEdd) l0g(?^Edd) 

Figure 14. Histograms of bolometric correction and Eddington ratio distributions. The grey filled histograms are obtained from the statistical method of host 
galaxy correction in the IR, whereas the hatch-filled histograms are from host galaxy SED fitting to remove IR galaxy contamination. 



5.1 Anisotropy and geometry corrections: smooth- vs-clumpy 
tori 

The similar regions of the Li2^m — ^2-iokeV plot in Fig. IH 
occupied by both obscured and unobscured sources suggests that 
the 12/im luminosities are isotropic (assuming that the absorption 
responsible for producing the MIR emission is broadly correlated 
with the X-ray absorption; see ^ and ^3.2.31 ). The SED models 
and anisotropy corrections u sed here in calculating Lbo i are based 
ultimately on the models of iGranato & Danes3 (Il994h which as- 
sume a smooth dust distribution in the torus. If the toru s is clumpy 
however, as recent models and fits to the data suggest (Ijaffe et aP 
l2004lTristrametal . 2007, Beckert et al. 2008), the torus emission 
is more isotropic and would not require significant anisotropy 
corrections between edge-on (obscured) and face-on (unobscured) 
AGN; this is possibly supported by the emergence of a unified 



— i^2-iokeV corrclation for both types of AGN and its im- 
plication of isotropic MIR emission. If we also inspect the values 
of 25//m luminosity plotted against X-ray luminosity we again 
find that obscured and unobscured objects do not occupy different 
regions of the plot, implying that the 25/im luminosity may also 
be, to a large degree, isotropic. In any case, even if we remove the 
isotropy corrections for obscured objects, this would have the effect 
of reducing bolometric corrections and Eddington ratios further 
by around 20-30 per cent for that subsample (since /^2-iokeV = 

(-t^X, total H~ /geo /anisotropy /starlight -t^ij^ ' ) /-t^2 — lOkeV r> 

for geometry, anisotropy and starlight correction factors 
/geo, /anisotropy , /starlight discusscd prcviously). Such a modifica- 
tion would not significantly alter our conclusions. 
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Figure 15. SEDs for the four objects with high bolometric corrections and Eddington ratios, under method 2 ( ^3.2.2^ for host galaxy removal. Key as in Fig. [8] 



5.2 Objects not detected by IRAS 

The non-detection of a substantial fraction of the sample warrants 
consideration of the expected distribution of values for A€2-iokeV 
and AEdd from the non-detected objects. Assuming the missing ob- 
jects are just below detection, the upper limits presented in Fig. [3] 
imply that the majority of the missing objects would have a very 
similar distribution of Li2/xm values to the detected objects. If the 
upper limits lay significantly above the relation in comparison to 
the detected sample, one would expect to see higher bolometric 
corrections on average due to a bigger potential nuclear IR contri- 
bution, and vice versa if they lay significantly below the line. Given 
that the upper limits also occupy a similar range in L2-iokeV as the 
detected sample, the distribution of bolometric corrections would 
also be expected to be similar. The mean black hole mass for the 
non-detected sample is log(MBH/M0) = 7.82 compared to 7.76 
for the detected sample, with very similar spreads, implying that the 
non-detections may lie at similar Eddington ratios as well. However 
this is obviously conjecture, and future work using IR data probing 
to lower fluxes would be able to verify or challenge these predic- 
tions. 



6 SUMMARY AND CONCLUSIONS 

We have presented a simple extension of the IPozzi et aP (l2007h 
method for calculating the bolometric output of AGN using repro- 
cessed IR emission from IRAS along with the hard X-ray emis- 
sion from Swift-ENI. We estimate the integrated nucle ar IR lumi- 
nosity by fitting of AGN nuclear SED templates from ISilva et al.l 
(I2OO4I) to the IR data, and estimate non-nuclear contamination in 
the large aperture IRAS fluxes via two methods: the first involves 
applying a correction to make the distribution of 12 /im fluxes lie 
on the most recent determination of the Lx — Li2^m relation de- 
termined bv lGandhi et al.l (l2009l) . and the second involves fitting a 
host galaxy SED template alongside the nuclear SED template. The 




41 42 43 44 

log(L2.iokev) from fit to BAT data 



Figure 16. Upper limiting infrared 12 iim luminosity against 2-lOkeV lu- 
minosity from fits to BAT data, for objects without IRAS PSC or FSC data. 
The completeness limit of the IRAS Faint Source Catalogue, 0.2 Jansky has 
been used in calculation of these upper limits. 

second method has the advantage that it does not presuppose a re- 
lation between X-ray and MIR luminosity, which will tend to bias 
the bolometric corrections to certain values, whereas the former 
method is superior in that it ought to give clean estimates of the 
core powers alone, on average. We fit a simple power-law model 
to the hard X-ray BAT data to determine the X-ray luminosities, 
and then calculate the bolometric luminosity after applying correc- 
tions for torus geo metry and anisotropy of emission as discussed by 
IPozzi et aP (l2007h . We calculate Eddington ratios using black hole 
m asses estimated f r om th e K-band bulge luminosity as described 
in IVasudevan et al.l (l2009l) (and use a well-constrained dynamical 



© 0000 RAS, MNRAS 000, 000-000 



18 R. V. Vasudevan, A. C. Fabian, P. Gandhi, Winter, L. M. & Mushotzky, R. F, 



mass for the known cD galaxy Cyg A), and present the 2-10 keV 
bolometric corrections. We calibrate our approach against the si - 
multaneous optical-to-X-ray SEDs from IVasudevan et"aD (l2009h 
and find trends for agreement between the two approaches, al- 
though there are still numerous sources of scatter. 

Using our zeroth-order IR and hard X-ray approach, bolomet- 
ric corrections for the low-absoprtion subsample give an average 
of around ^22, intermediate between the averages found from the 
two approaches of removing the host galaxy contamination (19,24). 
Eddington ratios are ^ 0.037 on average, from the two methods. 
For the high-absorption subsample the average Eddington ratios are 
0.034 and 0.043 from the statistical and host-SED fitting methods 
for host galaxy removal, respectively. The spread in bolometric cor- 
rections is greater than for low absorption objects, but for the lower 
part of the Eddington ratio distribution (AEdd < 0.03) bolometric 
corrections are consistently low in both methods for host galaxy 
contamination removal 15). The errors on these average bolo- 
metric corrections are typically less than ± ^ 3. Using the SED 
fitting method for host galaxy removal, for 0.03 < AEdd < 0.2 the 
high-absorption objects have bolometric corrections of 32 zb 10, 
possibly representing part of the transitional region from low to 
high bolom etric correct ions near Eddington ratios of 0.1 pro- 
posed by Va sudevan & F abian ( 2007), (2009), but better constraints 
on the host galaxy contamination and black hole masses are needed 
before this can be established. These bolometric corrections for our 
local, high-absorption AGN represent a new angle on uncovering 
the accretion processes at work in this important class of object. 
The objects not detected by IRAS {r^ 30 per cent of the potential 
sample) are expected to show similar results based on considera- 
tions of their upper limiting IR luminosities, but more sensitive IR 
observations are needed to confirm this. 

These preliminary suggestions that the majority of both high- 
and low-absorption //?A»S-detected objects in the Swift-EP^ cata- 
logue are at low Eddington ratios and exhibit low bolometric cor- 
rections reinforces unified scenarios in which the processes at work 
in AGN of different absorption levels are fundamentally similar, but 
the observed properties vary chiefly based on orientation. The pos- 
sible discovery of the lower end of an Eddington-ratio dependent 
bolometric correction for high- absorption objects would be an in- 
teresting aspect of this finding, subject to the uncertainties already 
discussed. 

Recently, iLamastra et al ] (l2009h have studied the Eddington 
ratio distribution of Seyfert 2s from the optical perspective, by esti- 
mating the optical bolometric correction factor to the [OIII] line lu- 
minosity. Fortheir sample of higher redshift AGN (0.3 < z < 0.8), 
they find that Seyfert 2s do not accrete close to the Eddington limit 
and Eddington ratios are ^ 0.1 on average. This work reinforces 
previous findings of predo minantly AEdd < 0. 1 accretion in the 
low-redshift universe (e.g. ISun&Malkai]ll989l) : here we extend 
this to an absorption-unbiased sample. 

High absorption objects are the dominant component of the 
X-ray background, and the identification of predominantly low 
bolometric corrections for them, using a larger sample than in 
IPozzi et aP (|2007) tells us how accretion onto black holes in the 
past can account for the c urrent popula tion of dormant massive 
black holes (as outlined bv lSoltanll 19821 and later authors). Lower 
bolometric corrections for low accretion rates found here are ap- 
propriate for matching the X-ray background energy density to the 
local black hole mass density ( Fabian 2004) . 

Additionally, the distribution of bolometric corrections in 
IPozzi et aP (l2007h (35 =b 9) is similar to the one we obtain for 
0.03 < AEdd < 0.2 sources (using method 2 for host galaxy cor- 



rection); indeed the Eddington ratios in the Pozzi et al. sample also 
lie predominantly in this range. Our extension of their method down 
to Eddington ratios as low as ^ 6 x 10~^ for a representative sam- 
ple of local obscured Seyferts is therefore valuable. 

This study is an extension of the Pozzi et al. technique of using 
the IR and hard X-rays (14-195 keV) for estimating the bolomet- 
ric output in a way that is less prone to the spectral complexities 
which affect the optical-UV and X-ray regimes, particularly for 
obscured objects. However, as discussed at length, appropriate cor- 
rection for the host galaxy contamination in the IR is essential. A 
useful extension of this study would be to use higher resolution 
Spitzer data in combination with BAT data to perform the same 
calculation, which would provide a more accurate measure of the 
nuclear reprocessed IR without hav ing to apply the large contami- 
nation corrections needed here fe.g. lMelendez et al.l2008h . The un- 
certainties in the covering fraction also represent a limitation of the 
study, and future work determining the appropriate distribution of 
covering fractions to use will be very valuable. 

Some of the uncertainties in the black hole mass estimates 
should be addressed in the analyses of Winter et al. and Koss et 
al. (in prep) which will provide Hp linewidth-based estimates of 
the black hole masses in the catalogue, which can be used to refine 
the position of the sources on the Eddington ratio axis. 

One class of object for which our approach may not be ap- 
propriate is the Ultraluminous Infrared Galaxy (ULIRG) class. In 

these, the mid-infrared continuum is not directly attributable to the 
ji 1 

Alonso-Herrero et al . 2006) and so they should be cor- 
rected for or excluded when consid ering the IR continuum of large 
samples of AGN, as discussed by [Thompson et al.l (l2009h . How- 
ever the redshift distribution of ULIRGs significantly differs from 
the BAT catalogue, so we expect contamination to be minimal, 
with none of the objects presented in our study being classified as 
ULIRGs in NED. 

Questions remain as to exactly what the IR continuum due 
to reprocessing in AGN is. The continuum produced depends on 
the dust distribution, grain size and torus geometry; in particular 
the dumpiness or smoothness of the dust distribution is impor- 
tant. Another useful extension to this work would be to explore 
how different continuum models, for example those for clumpy tor i 
dHonig et al.ll2006L iNenkova et al.ll2008aL iNenkova et al.ll2008bb . 
would alter the values of Lboi found here, especially as recent ob- 
servational evidenc e strongly favours the clump y torus scenario. 
The recent work of iRamos Almeida et al ] (l2009h . who fit clumpy 
torus models to mid-IR SEDs for a small sample of AGN, paves 
the way for further studies in this area. 

The 9-month BAT catalogue dTueller et al.l2008h has provided 
numerous opportunities to understand the accretion processes at 
work in local AGN (e.g. iMushotzkv etaDl2008L IWinter et al.ll2009l . 
IVasudevan et al ]|200i). Similar follow-up studies on the recently- 
released 22-month BAT catalogue will afford the chance to widen 
these analyses to a larger, more representative sample of local 
AGN. 
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